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ABSTRACT 

The differential equations governing the propagation of sound in a 
variable area duct or nozzle carrying a one-dimensional subsonic 
compressible fluid flow are derived and put in state variable form using 
acoustic pressure and particle velocity as the state variables. The duct or 
nozzle is divided into a number of regions. The region size is selected so 
that in each region the Mach number can be assumed constant and the area 
variation can be approximated by an exponential area variation. 

Consequently, the state variable equation in each region has constant 
coefficients. The transmission matrix for each region is obtained by 
solving the constant coefficient acoustic state variaole differential 
equation. The transmission matrix for the duct or nozzle is the product of 
the individual transmission matrices of each region. Solutions are 
presented for several geometries with and without mean flow. 

NOMENCLATAURE 



matrix 

A 

area 

B,C,D 

matrix 

b,c 

constants 

c 

isentropic speed of sound 

f 

frequency, Hz 


I 


2 


I unit matrix 

i 

t<Q propagation wave number, ui/Cq 

Mq Mach number, 

(MW) molecular weight 

m d ln(A)/dx 


P transformation matrix , , 

p pressure % 

R gas constant 

s entropy of gas 

T transmission matrix 

t time 

u fluid velocity 

M mass flow , 


X Cartesian coordinate 

Y acoustic state vector 

Z acoustic impedance 

AX region length 

Y specific heat ratio of gas 

e m/2ik^ 

f Ko/(l - 

X . eigenvalue 

e temperature 

Pq gas density 

u) angular frequency, radians/sec 

Superscripts and subscripts: 

( ”* ) vector quantity 
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( ~ ) instantaneous quantity 

( )q reference state quantity 

( perturbed quantity 

INTRODUCTION 

The use of the transmission matrix approach provides considerable 
simplification to the algebraic complexities involved in the analysis of 
acoustic systems which can be modeled as a series of simple regions. The 
elements of a region transmission matrix are a set of coefficients which 
depend only on the properties of the region. Furthermore, this set of 
coefficients relate the acoustic pressure and particle velocity at the inlet 
of the region to those quantities at the exit of the region. The algebraic 
complexities are reduced in this case because the transmission matrix of the 
entire system is the product of the constituent transmission matrices. 
Consequently, the properties of a system can be calculated by 
straightforward multiplication using a building block approach based on 
defining the input-output relationship of each of the regions of a system. 

nV' '■ 

As part of a combustion noise research program conducted at the NASA 
Lewis Research Center, internal pressure measurements are made in aircraft 
engines. Because the frequency spectrum of the combustion noise peaks below 
1000 Hz and the combustion noise generally propagates in a small area 
annulus, the plane wave mode of propagation is extremely important. Plane 
wave combustion noise propagation is a ducted combustion system with flow 
has been studied by Miles and Raftopoulos (refs. 1 and 2). Miles and 
Raftopoulos constructed a system model using a transmission matrix approach 
to the relate acoustic pressure and particle velocity at one point in the 
system to those quantities at another pointii^ In constructing the model, it 
was assumed that the duct system consisted of constant area duct regions 
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which could be matched at area discontinuities. Then transmission matrices 
were calculated for both the constant area duct regions and the area 
discontinuities using the systems acoustic differential equations. 

This approach of using transformation matrices should also be 
applicable to the study of the propagation of plane wave combustion noise in 
aircraft engine ducts and nozzles. However, aircraft engines do not have 
only constant area ducts. Also, by definition nozzles have varying area. 
Thus, as a first step in constructing an aircraft engine nozzle or duct 
system model for combustion noise propagation, this paper derives the 
input-output relationship in the form of a transmission matrix for the 
propagation of sound in a variable area duct or nozzle carrying a 
one-dimensional subsonic compressible fluid flow. 

BACKGROUND 

Many studies of acoustic problems related to the acoustics of aircraft 
engine-duct systems have been made. Much of the previous work was related 
to fan noise and this type of work>is reviewed in reference 3. While the 
acoustics of variable area ducts is among the topics reviewed in 
reference 3, the transmission matrix approach is not discussed in 
reference 3. This topic was probably omitted because the transmission 
matrix approach is used to study the plane wave propagation of low frequency 
noise and is not used for fan noise. However, the application of 
transmission matrices to ducts without flow is discussed in references 4 
to 11. Among these references, note that a brief but comprehensive 
introduction to the fundamentals of the transmission matrix representation 
is given by Lampton (ref. 9). As mentioned previously, transmission matrices, 
can be used to model duct systems with flow if the systems consist of 
different diameter constant area ducts by taking into account the area 
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di 5Contiiiuities. Methods for dealing witlj area discontinuities are 
discussed in references 1 and 12 to lb. 

To the author's Knowledge, transmission matrices have not been used to 
model systems containing variable area regions with mean flow. However, if 
there is no mean flow in a variable area duct, the acoustic equations reduce 
to the Webster Horn equation for which solutions may oe obtained for some 
specific duct shapes (refs. 27 to 21). Consequently, it is possible to 
determine exact transmission matrices for these cases. 

With a mean flow in a variable area duct, the problem of finding a 
transmission matrix becomes more complex, since the velocity, density, and 
Mach number are then dependent on the duct area. Consequently, the acoustic 
differential equations have variable coefficients. An exact solution to the 
differential equations governing the propagation of sound in a variable area 
duct with one-dimensional flow was obtained by Eisenoerg and Kao (ref, 22). 
However, this solution is hot applicable to the present problem since it was 
for specific duct shapes which are semi-infinite and diverge from a sonic 
throat. Furthennore, the equations are not in transmission matrix form. 

For ducts having more general shapes, numerical methods can be used to solve 
the variable-coefficient equations. Solutions of the one-dimensional 
acoustic differential equations with variable coefficients have been 
obtained by Davis and Johnson (ref. 23), King and Karamcheti (ref. 24), and 
Lumsdaine and Ragab (ref. 25) using numerical methods. Again none of these 
solutions is in transmission matrix form. 

Another numcrn'cal technique for calculating the behavior of sound is an 
approximate technique developed by Alfredson (ref. 26) which uses many 
small discontinuities to define the duct shape. Each discontinuity is then 
treated using the procedure described in reference 12. The technique 
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described herein also uses many small discontinuities to describe the duct 
shape. However, trie method is based on the infinite exponential horn 
solution (described in refs. 17 and IB) and includes mean flow effects. 

The governing equations used herein have been used by H. S. Tsien 
(ref. 27), L. Crocco and S. Cheng (ref. 28) to study oscillations in rocket 
nozzles, and S. Candel (ref, 29) to study entropy noise produced by a nozzle. 

In this paper, first, the governing equations are derived by 
linearization and, put into state variable form using acoustic pressure and 
particle velocity as the state variables. Next, ih order to solve the 
resulting differential equation, tne duct or nozzle is divided into regions 
so that the coefficients of the differentia! equation can be made constant. 
The transmission matrix is obtained for each region by solving the acoustic 
state variable differential equation using methods oescribed by Ogata 
(ref. 30). The transmission matrix for the variable area duct or nozzle is 
the product of the transmission matrices for each segment. Last, solutions 
are presented for several geometries with and without mean flow. 

DERIVATION OF STATEI VARIABLE EQUATION 
The following assumptions are inade* The duct or nozzle has an area 
profile A(x). The flow in the duct or nozzle of an ideal nonviscous, 
nonconducting perfect gas is described by the following equations of 


continuity, momentum, energy, and state: 

a(pA)/3t + 3(puA)/3X = 0 (i) 

p(3u/3t + usuT/ax) + 3P/3X = 0 ■ (2) 

(3/3t + ua/ 9 x)s =0 (3) 

S - Sq = C^ ln(p/p'^) (4) 


Where p, p“, u, and s' are, respectivelyj the instantaneous pressure, 
density, velocity, and entropy. 


The fluctuations in the duct dr« assumed so small tnat the flow in the 
duct 1s only slightly perturbed. Consequently^ the instantaneous quantities 
can be written in terms of a unperturbed stationary flow quantity designated 


by 0 and a small perturbed quantity designated by i as follows: 


P ** Pq(x) * Pj^(t,x) (5 

e « Po(x) + Pi(t,x) (6 

u *= Uq{x) + U3^(t,x) (7; 

$ a Sq(x) + s^(t,x) (8; 

Substituting equations (6) to (8) into equations (1) to (4) yields the 
following zeroth order system of equations: 

d(PoUQA)/dx 0 (9; 

PqUq du^dx + dp^/dx s 0 (10] 

Sq s Cy Inip^/pJ) - constant (11) 

and the first order system of equations: 


A 3p^/3t + u^A 3pj/3X + APq 3u^/3X + p^(Uq dA/dx + A du^/dx) 


+ u^(p^ dA/dx + A dp^/dx)= 0 (12) 

pQ 32/3t + PqU^ QU^/dx + PqU^ auj^/3x + p^u^ du^/dx + 3p.j^/3X = 0 (13) 

(3/3t + 3/3X)s^/C^ = (3/3t + U^ 3/3X)(p^/p^ - YPj/Pq) « 0 (14) 

For the system under consideration the gas equation is 

p = F^?o/(MW) (15) 

and the isentropic speed of sound is 

. t«8o/(MW) =, rp„/p„ (16) 

From equation (14) the density perturbation is related to a pressure 
perturbation by 


Substituting equation (17) Into aqutitlons (IH) and (13) ylelas the 
follov/ing equations 1n matrix form: 


0|/apj/at\ ll/»Pl/9x\ 




Taking the Fourier transform of equation (18) and solving for dY/dx 


yields 


dY/dx = BY 


Using the zerotn order system of equations (y) to (11) and equation (16) to 
relate the gradients of the logarithmical variation of velocity and density 
to the gradient of the logaritnniical variation of area yields 


Using equations (20) and (21) the elements of the B matrix in 
equation (19) are 



y 


B 


M 


* (1 - M^) 


0 ti 1n(A) . 


(I - %) 


12 


^ pQ^o%^ (I in(A) 


«\2 CIX 


ik. 


oi ® 


11 


- My) 


(22) 

(23) 

(24) 



<1 " >^o) a ln(A) . "ot-^y 

(l - »lf (1 - 


The acoustic state variable differential equation is given by 


equation (19). This equation will now Oe solved to obtain the duct 
transmission matrix, 

METHOD OF SOLUTION 


(25) 


In the previous section a general acoustic state variable differential 
equation was derived (eq. (19)). In this section an approximate solution to 
equation (19) is obtained. The first step in solving equation (19) is to 
divide the duct or nozzle into a number of regions or subsections. The 
region size is selecteo so that the Mach number can be assumed constant in 
each region and the area variation in each region can be approximated by an 
exponential area variation. Consequently, in any given region 

Mq =5 constant (26) 

A ^ A^e'"^ (27) 

w'ith these assumptions the Mach number factors in the B matrix in 
this region are constant. Also, the value of d(1nA)/dx is given by 
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0(1 n(A)) 
ux 


.1 HI 


(;^B) 


The B matrix is now independent of x in the region where it is to 
be evaluated. Thus, the solution to equation (19) for the region is 


Y(Xj + AX) « {e'^^^)j.Y(Xj.) 

ith 


(29) 


where ax is the length of the region and exp(BAx) is known as the 
matrix exponential. The value of exp(BAx) is defined by 


Jax 

e ss 


I + (Bax) -i- + . . . 


IT 


(30a) 


The matrix exp(BAx) is the transmission matrix of the region, 


(T)j, so that 


, (Bax). 
(T)j = e J 


00 

^^(Bax)^ 


f=0 




'21 


12 


T. 


22 


(30 b) 


Consequently, the transmission matrix, T, for the variable area duct is 
found from 


'll 


Y(x) 


^TT(e'^^)^Y{x) 

x=L I I J 


j=l 


X=;0 


T. 


21 


^121 


”'’221 


Y(X) 


X=0 


= TY(x) 


lx=0 


(31) 


where 


L = 2 (AX) 


j = l 


(32) 
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and N Is the number of duct regions. 

In this section an approximate solution to the state differential 
equation was obtained, In the next section computational techniques useful 
for obtaining analytical and numerical solutions will be discussed, 

EVALUnON OF mTRIX EXPONENTIAL 


The solution of the matrix space differential equation thus involves 
the matrix exponential exp(Bax), It is necessary to express exp(Bax) as a 
matrix in order to obtain numerical values for the transmission matrix. 
Methods for evaluating a matrix exponential are discussed by Ogata in 
chapter 6 of reference 30, These same methods can be used for evaluating 
exp(Bax). The first method discussed produces a solution which can be 
evaluated on a digital computer using the built-in complex exponential 
function. It is also useful in generating analytical solutions in terms of 


exponential and trigonometric functions. The second method that will be 
discussed uses equation (30) to numerically obtain a value of (T) by 
evaluating the exponential power series in (Bax). 

The first methoo is to transform the B matrix into a diagonal matrix 
D using a transformation matrix P. If the eigenvectors of the B matrix 
are distinct, then 


where 









(33) 



Tti!s method Is implemented by first calculating 


so that 


si m CBC 


-1 


'"^12 ** 

‘^21 “ 
sir)!) ** 

W MM 

Next, the eigenvalues of the si matrix, A, are found by setting the 
determinant of the matrix (j/ - aI) equal to zero. Consequently, tne 
eigenvalues are the roots of the equation 


A 



bA + c = 0 


Where 

b « '*■ -*^22^ 

and 

c = j/il-^22 " •^12'^21 

Next, the P transformation matrix is calculated using 
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P « 


"^11 
•*^11 “ h 


“•*^21 

j/ll - 


(41) 


Note that if the P matrix defined by either equation (40) or (41) is 
multiplied by a constant the result is also a P matrix. Consequently, it 
is possible to significantly simplify the P matrix generated by equation 
(40) or (44). The diagonal D matrix is then 







Consequently, the square matrix, exp(Dax) is given by 


(42) 



eOAX 


Xj^AX 

e 


0 


0 

XoAX 
e *■ 


The value of the square matrix exp(Bax) is then given by 


(43) 


11 


12 1 


T = e^^^ = C"^Pe*^^^P~^C = 


(44) 


1 21 


*22 


Note that the elements of the matrix contains only complex exponential 
functions from equation (43). 

As mentioned previously, the second method for evaluating exp(BAx) is 
based on expanding exp(BAx) in the exponential power series in (Bax) given 
by equation (30). Equation (30) is also important in obtaining analytical 


solutions. Fon instance, using equation (30) it is possible to prove 
equation (43). The virtue of equation (30) is its simplicity. However, in 
writing a digital computer program based on equation (30) convergence 
requirements must be specified and this may increase the running time of a 
computer program evaluating exp(Bax). 

In this section two procedures for implementing the solution were 
presented. Some analytical solutions will be presented next. These 
solutions are presented to clarify the computational techniques and to 
provide some insight into the numerical results. 


ANALYTICAL SOLUTIONS 

In this section analytical solutions for the input-output relationship 
of the acoustic pressure and the particle velocity in the form of a 
transmission matrix will be obtained using the methods described in the 
previous section. The solutions are exact since only cases where the Mach 
number is zero or constant are considered. The geometries considered are a 
duct with constant area and a duct with an exponential area variation. 

These solutions are presented to clarify the numerical results presented in 
the next section and the methods used to obtain them. 

^ ■ ' l:/ ' 

Constant Area Duct Without Flow 

First, consider a constant area duct of length L with no flow. For 
this case M = 0 and d ln(A)/dx = 0. Cons^equently, substituting these 
values of M and d 1n(A)/ax into equations (22) to (25) yields the 
following B matrix 



0 



ik^/p„c„ 

0 0 0 


p„c„ik„ 

0 0 0 




Substituting equations (45) and (32) into equation (31) yields 
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since the B matrix is independent of the Mach number. Using equation (30) 
to evaluate equation (46) yields 



(47) 


Substituting the series definitions of the trigonometric functions, it 
follows that 



cos(k|jL) 


i; sin(l<QL) 


sin(k„L) 
0 0 'o' 

CQs(k^L) 



Constant Area Duct With Mean Flow 
Next, consider a constant area duct of length L with mean flow. 

Since the duct area is constant again d ln(A)/dx = 0. Using equations (22) 
to (25) the B matrix is given by 



1 

! 


- 

0 

1 

* "o 


M,(-ik,) 



1 - 

"o . 


Substituting equation (49) into equation (31) shows that the transmission 
matrix is again given by equation (46). For this case the value for 
exp(BL) will be evaluated using equation (33). 

First, the matrix is calculated from equation (35). .Next, the 
eigenvalues are found using equations (37) to (39). Thus, 

*r I 

*2 * " *^0* J 

Based on equation (40) a satisfactory P matrix is 

n -n 


The diagonal D matrix calculated by substituting equations (34), (35) 
(40), and (51) into equation (42) is | 


The transmission matrix calculated by substituting equations (34), (51), 
and (52) into equation (44) is 


X-iL XnL] 


x^L XgLj 
.e - e / 


X,L XpLl 
le - e ) 


X,L XnL] 

.e + e J 


or using the complex definition of the sine and cosine functions 


:co${?|.) ip^Cg sin(tL) 


T ^ e 


-itwa 


i sin(gl] 


cos(tL) 


where t = kj{l - Mq). 

Note that for M = 0, the transmission matrix given by equation (54) 
reduces to that given by equation (48), 

Variable Area Duct Without Flow 

The second geometry discussed in this section is a duct of length L 
having an exponential area variation 


A - A e"*- 
Ml - A^e 

First, tne zero frequency limit transmission matrix for no flow will be 
calculated. For this case M = 0, k = 0, and d ln(A)/dx = m. 
Consequently, from equations (22) to (25) the B matrix is given by 
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(b6) 


0 -nt 


Substituting equation (86) into equation {31) again yields equation (46). 
For this case, equation (46) can be evaluated using equation (30). Again 
substituting for tiie power series obtained yields a closed form solution as 
follows 


T = 


1 


0 


,-mL 


Substitution of equation (55) into equation (57) yields 


T = 


0 A„/Aj 


The corresponding set of simultaneous equation is then 


(87) 


(58) 


p(ia),L) = p(io),o) 


(59) 


A|_u(iio,L) = AyU(ia),o) (60) 

Consequently, at low frequencies for no flow the acoustic pressure is the 
same at both ends of the duct and the acoustic volume flow is constaht. 

Next, the general transmission matrix for the exponential duct without 
flow will be calculated. For this case d ln(A)/dx = m. Consequently, from 
equations (22) to (25) the B matrix is given by 


Substituting equation (61) into equation (31) snows that the transmission 
matrix is given by equation (46). For tnis case, the value of exp(BL) will 
be evaluated using equation (33). The eigenvalues found using equations 
(37) to (39) are 



The diagonal D matrix calculated by substituting equation (64) into 
equation (42) and using equation (35) to determine the sf matrix is 



The transmission matrix calculated by suostituting equations (34), (63) 
and (65) into equation (44) is 






n 


Note that for a constant area duct ni » 0 and thus e 0. The 
transmission matrix given by equation (67) for this case reduces to the 
transmission tnatrix given by equation (46). 

Also note that for the low frequency case, the eigenvalues are given by 

lim » 0 I 

k»0 ^ 

> ( 68 ) 

11m Xo « -m I 

k>0 J 

Substituting these eignevalues into the transmission matrix given by 
equation (66) and letting k go to zero, yields the transmission matrix 
given by equation (57) . 

In this section analytical solutions to equation (3i) were found for 
the case where the Mach number in equations (22) to (2b) is zero or constant 
and equation (31) yields an exact solution. It was shown that using 
equation (33) to obtain a solution for exp(BL) yields a solution in terms 
of exponential or trigonometric functions. It was also demonstrated that 
while equation (30) can be used to obtain an analytical solution. It is only 
useful in the most simple cases where the resulting power series 
recognizaule. The next section discusses numerical solutions. 

COMPUTATIONAL EXAMPLES 

Numerical calculations of transmission matrices for several duct 
configurations with and without flow were made by evaluating equation (31). 
The calculations yield exact solutions for cases where the Mach number is 
zero or constant. Tnus for any of the cases discussed above the numerical 
solutions are identical to the corresponding analytical solution. 

When the Mach number is not constant or zero the numerical calculations 
yield an approximate solution. The approximation is improved as the duct is 
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divided into smaller sudsections. Typical results obtained by dividing the 
duct configuration into 500 equal length subsections are presented next. 

For all solutions, the zeroth order system of equations is solved first to 
evaluate pg, Uq, ©o* ^o the start of each subsection using a 
stagnation pressure of one atmosphere and a stagnation temperature of 
295 K. The ratio of gas specific heats was 1.4. 

The numerical results are presented as plots of the magnitude and phase 
angle of each of the matrix elements versus frequency. 

Constant Area Duct 

Results obtained for a one-meter- long cylindrical duct will be 
discussed first. Since in a cylindrical duct the Mach number is zero or 
constant, the numerical solution is identical to the analytical solution. 
Transmission matrices for the one-meter-1 ong constant area duct calculated 
for M - 0, M = 0.25, and M = 0.6 are shown in figure 1. The solid curves 
in figure 1 correspond to solutions obtained using equation (48). The 
dashed curves correspond to solutions obtained using equation (54). The 
curves are functions of = (2Trf/c^)L/(l - M^). Consequently, the 
curves appear to contract with increasing Mach number. This is because the 
frequency location of corresponding peaks and dips which occur at constant 
values must shift to lower frequencies as the Mach number is increased. 
In addition, as the Mach number increases so does the phase factor 
exp(-i£MQL) which multiplies each transmission matrix element. 

The ratio of the pressure at x = 1.0 to the pressure at x = 0 with 
zero acoustic particle velocity at x « 0, T^, is equal to 
exp(-ifMjjL)cos(f L) which is plotted in figure 1(a). The ratio of the 
acoustic particle velocity at x = 1.0 to the acoustic particle velocity 
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at X a 0 with i^ero acoustic pressure at x » 0, T^ 2 » shown in figure 
1(d) and is also equal to exp(-itMpL)cos(^L) . 

The ratio of the pressure at x « 1.0 to the acoustic particle 
velocity at x « 0 with zero acoustic pressure at x « 0, T^ 2 » 

^ 0*^0 - &MQL))sin(fcL) which is shown in figure 1(b). The 

ratio of the acoustic particle velocity at x = 1.0 to the acoustic 
pressure at x = 0 with zero acoustic particle velocity at x « 0, T 2 ]^, 
is exp(i(ff/2 -- tMQL)sin(tL)/pQCg which is shown in figure 1(c). 

Variable Area Ducts 

Results obtained for one-meter- long ducts having exponential area 
variation will be discussed next. Two examples are presented. One example 
which is identified as the long exponential nozzle (LEN) has an area which 
decreases exponentially from 0.02 to 0.01 m (d ln(A)/dx - = ln(l/2)). 

The transmission matrix for this geometry calculated for mass flows of 0 and 
2.0 kg/sec are shown in figure 2. The second example which is identified as 
the long exponential horn (LEH) has an area which increases exponentially 
from 0.01 to 0.02 m (d ln(A)/dx ^ m^^j^ = ln(2)). The transmission matrix 
for this geometry calculated for mass flows of 0 and 2.0 kg/sec is shown in 
f i gure 3 , 

Without flow . - First, the no-flow cases will be discussed. At low 
frequencies, the analytical solution given by equation (58) for the long 
exponential nozzle transmission matrix, ^LEN» 


- 0 ) = 


[o (%/^l)lenJ 


(69) 
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whefe exp(-in^gj^L) * * 2-0. Furfchennore* at low 

frequencies, the analytical solution given by equation (58) for the long 

a 

exponential horn transmission matrix, T, is 


3- 


\eh 


1 


0 


(70) 


where exp(-iii|_E„L) . (A„/A|_>le„ - 1/2. 

In addition, at frequencies greater than |mcQ/4w| the analytical 
solution given by equation (67) for is 


'LEN 



1/2 


LEN 


cos(KqL) IpqCq sin(kQL)’] 


i sin(k^L) 

'"p c ” 

*^0 0 


cos(kQL) 


(71) 


where (Aq/Al)l£Hj = 2. Moreover, at frequencies greater than 
jmcQ/4irj the analytical solution given by equation (67) for is 


T. 


LEH 



1/2 


LEH 


cos(k^L) 


i sin(kQL) 


iPoCo sin(k^L) 


P^nC^ 
0 0 


cos(k^L) 


(72) 


where (Aq/Al)^^^ = 1/2. The results of the numerical calculations 

shown in figures 2 and 3, for the case without flow are identical with these 


analytic results, These few exatnples show that the numerical calculations 
are in agreement with the analytical solutions, 

With flow , - The numerically calculated transmission matrix for the 
long exponential nozzle and horn with a mass flow of 2,0 kg/sec shown in 
figures 2 and 3 will be discussed next. 

As expected from the results presented in figure 1 the frequency 
location of the peaks and dips shifted to lower frequencies. Also as 
expected, the phase angle is significantly different for this case. 

However, the types of changes in the magnitude and phase angle of the 
elements of the long exponential nozzle and horn transmission matrix is 
unexpected. The interaction of the flow and the gradient of the 
logarithmical area variation change the transition matrix elements 
significantly. For Instance, the Vow frequency value of {T 22 )Lgf^ 
shifts from 2 with no flow to 3.3 with a mass flow of 2 kg/sec. Also, the 
low frequency value of shifts from 1/2 with no flow to 0.2y5 

with a flow of 2.0 kg/sec. 

The numerical method can be applied to any type of area variation. The 
transmission matrix for a short exponentially contracting duct and a 
linearally contracting duct, both 20 cm long, are compared in figure 4. 

Below 50 Hz the transmission matrices are similar. Above 60 Hz the 
magnitudes do not agree. However, the phase angle differences become 
significant only above 400 Hz. 

CONCLUDING REMARKS 

A numerical method for detenaining the transmission matrix of a 
variable area duct or nozzle carrying a compressible subsonic mean flow was 
presented, The method yields an exact solution for the zero mean flow case 
and an approximate solution for the nonzero mean flow case. However, in 
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using the numerical method the duct or nozzle Is divided Into regions in a 
way that makes the method more accurate as the number of regions Increases 
and the region size decreases. Consequently, the accuracy of the method 1s 
only limited by the region size. 

Furthermore, for the no-flow case not only tloes the numerical solution 
become exact, but It also becomes possible to use the equations which 
describe it to generate an analytical solution in terms of exponential and 
trigonometric functions. 

The method was applied to several geometries with and without flow, 
1n<”luding both expanding and contracting duct sections. For the cases 
without flow analytical solutions were obtained. 

The numerical method may also be useful In studying sound propagation 
in ducts without a mean flow having shapes for which the Webster Horn 
equation does not have an exact solution. Furthermore, the method is easy 
to apply. Therefore, even for ducts with shapes for which the Webster Horn 
equation does have an exact solution, the method may be useful where a 
numerical solution is adequate. 
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Figure 1. Cylindrical duct transmission matrix, T, 
where YU • 1. 01 • T?U ■ 0 OML • I m, T ■ 295 K, 
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Figure Z • Long exponential nozzle transmission 
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Figure 3. - Long expqneniial Hornjrinstnission 
melrix, T, where Y(x • 1. 0) • T?(x • 0. 01. 
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